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INTRODUCTION 


Rising energy costs as well as current and proposed noise regulations have 
prompted continuing research and development of lighter weight, improved noise- 
suppression systems for the turbofan engines of commercial aircraft. The use of 
uniform acoustic duct liners (liners with wall impedances that do not vary axially 
or circumferentially) to suppress the noise generated by these turbofan engines is 
well established. However, the above factors require that acoustic-liner designs be 
continually refined and updated to incorporate new concepts. Recently, variable- 
impedance liner concepts (fig. 1) have been considered in the design of more effec- 
tive noise-suppression systems. 

Literature on axially segmented liners, a concept initially derived by Zorumski 
(ref. 1), is exhaustive. Lester and Posey (ref. 2) have shown that an axially seg- 
mented liner with two sections gives more suppression than the best uniform liner. 
Sawdy, Beckemeyer, and Patterson (ref. 3) have shown that this increase in suppres- 
sion is due to a scattering of the sound at the first axial discontinuity into higher 
order radial modes, which are then more highly damped by the second section. How- 
ever, certain characteristics of axially segmented liners limit their usefulness. 

These include lack of a broadband attenuation spectrum and inability to outperform 
uniform liners for well cut-on modes (ref. 4) . 

Mani (ref. 5) observed experimentally that an additional suppression of as much 
as 5 dB could be obtained at some directivity angles by circumferentially segmenting 
a uniform liner. He suggested that the increased suppression was because of scatter- 
ing of acoustic energy into higher order circumferential modes, which were more highly 
damped. In addition, Mani suggested combining axially and circumferentially seg- 
mented liners to form a checkerboard pattern. These experimental findings motivated 
researchers to develop analytical models to verify these findings. Thus, in this 
instance, duct-acoustics experimental work has preceded and guided analytical 
development in liner technology. 

Development of reliable analytical models for predicting noise suppression in 
ducts with circumferentially segmented liners is a prerequisite to the optimum design 
of these liners. Watson (ref. 6) developed a no-flow model for a rectangular duct 
based on a finite-element method, while Astley et al. (ref. 7) incorporated the 
effects of a sheared mean flow into a circular duct using the finite-element method. 
Since the finite-element methods developed in references 6 and 7 are in core solu- 
tions, the number of circumferential strips is severely restricted and the effects of 
the sound source are difficult to evaluate. These models are therefore not suited 
for optimization studies. Namba and Fukushige (ref. 8) used the equivalent surface- 
source method to evaluate attenuation properties of circumferentially segmented 
liners in a rectangular duct with uniform mean flow. They concluded that circum- 
ferentially segmented acoustic liners bring about modal transfer of sound energy 
between modes of different circumferential wave numbers and lead to higher power 
attenuation levels and broader effective attenuation bandwidth than a uniform 
liner. Watson (ref. 9) used the Galerkin method (ref. 10) to evaluate attenuation 
properties in circular ducts without mean flow. Results indicated that circumfer- 
entially segmented liners scatter energy between circumferential modes of different 
orders and lead to broader attenuation curves. However, before circumferentially 
segmented liners can be fully evaluated, additional studies evaluating the effects of 
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frequency, source complexity, and duct length are needed. Such liners should also be 
optimally designed to achieve maximum suppression and their off-optimum properties 
should be explored before their potential benefit can be fully realized. 

This paper assesses the relative merits of circumferentially segmented liners 
by comparing their optimum attenuation properties with those of optimized uniform 
liners. Off-optimum properties are also discussed as well as effects of frequency, 
source complexity, and duct length. The analytical prediction program developed in 
reference 9 is used in conjunction with an optimization method to compute optimum 
properties for the uniform and segmented liner. The optimization method uses the 
Davidon-Fletcher-Powell algorithm (ref. 11). However, a series of two-dimensional 
contour maps is used to obtain good starting values for the algorithm. Both plane 
and circumferential standing-wave sources are considered in the study. Finally, 
emphasis is placed on segmented-lining configurations consisting of hard-wa 11/soft- 
wall combinations of the kind considered by Mani (ref. 5). 


ANALYSIS 

Consider a semi-infinite circular duct with radius a, as shown in figure 2. 

(A list of symbols used in this paper appears after the references.) A circumfer- 
entially segmented liner with specific acoustic admittance p(6) is placed along the 
outer wall, and an initial noise source in terms of a pressure boundary condition is 
given at z = 0. Further, the segmented liner is assumed to consist of two piece- 
wise uniform liners with admittances and ^ 2 ^ shown in figure 3. These 

piecewise uniform liners are combined in circumferential series to form a total of 
2T strips. The admittance function p(0) is expanded as the Fourier series 


^2^ • . sin 3T0 . sin 5T0 . 

= Pav ^ ^ ^ 


...) 


( 1 ) 


where is the average value of p(0) for the circumferentially segmented liner 
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p(0) de = 




and T is the periodicity of the liner. Within the context of this paper an opti- 
mized circumferentially segmented liner is defined as a liner configuration possess- 
ing uniform admittance p^ and p^ chosen to maximize the transmission loss over 
the length L. Finding this optimized liner is the major thrust of this paper. 

Dimensionless acoustic waves propagating within the duct depicted in figure 2 
satisfy the Helmholtz equation (ref. 9) 

P(r,e,z) + P(r,e,z) = 0 (2) 


in which a time dependence of the form e"^*^ has been assumed, K = wa/c is 
the wave number, P is the dimensionless acoustic pressure, co is the angular 
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frequency, V is the Laplace operator, t denotes time, and c is the ambient 
speed of sound. Throughout this paper, all distances are referred tg the duct 
radius a and the acoustic pressure is referred to the quantity pc t in which 
p is the ambient density of the medium. In addition, the boundary condition along 
the outer wall is expressed in terms of the admittance p(0) as follows: 

= iK p(0) P(1,9,z) ( 

whereas the boundary condition at the source plane is expressed as follows: 


P(r, 0,0) = Q(r, 0) 


(4) 


where Q(r, 0) is the known source pressure function. 

Equations (2) to (4) constitute a boundary value problem for the acoustic pres- 
sure field P(r,0,z). Furthermore, this field can only contain waves running to the 
right because the duct in figure 2 is semi -inf ini te. In this paper, the solution is 
obtained by the eigenfunction expansion technique (ref. 3). The eigenfunction expan- 
sion technique expands the acoustic pressure field in terms of the eigenfunctions of 
the circumferentially segmented liner: 


P(r,0,z) = ^ A P (r,0) ej^(iK z) 
Z=0 ^ ^ ^ 


(5) 


in which Kj^ is the complex axial wave number and P^(r,0) is the acoustic pressure 
eigenfunction for the circumferentially segmented liner. Solutions for the eigen- 
function axial wave number are developed in the following sec- 

tion. These eigenfunctions can also be shown to be orthogonal (ref. 9), so that 
substituting equation (5) into equation (4) and making use of the orthogonality con- 
dition allows the following determination of A^: 


A 


Z 



Q(r, 0) P^(r, 0) r dr d0 



P^(r,0) 


r dr d0 


( 6 ) 


A measure of the liner effectiveness as a sound absorber is given by the 
transmission-loss function TL as 


TL(p^,P 2 > = 10 log^^[W(0)/W(L)] 


(7) 
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and 


W(z) 



iP* 


5P 

dz 


r dr del 


( 8 ) 


where Re{ ) denotes the real part of the complex expression within the parentheses 
and the superscript asterisk indicates the complex conjugate, Ihe value of the 
transmission-loss function can be obtained once the solution for the eigenfunction 
and for the axial wave number is obtained. The methodology used to obtain these 
solutions is described in the following section. 


Solution for Eigenfunction and Axial Wave Niomber 


If equation (5) is substituted into equations (2) 
(r, 9) can be shown to satisfy a Helmholtz equation 


and (3), each eigenfunction 


^ Pj^(r,0) + Pj^(r,e) = 0 


(9) 


with the homogeneous boundary conditions 


5P^(r,0) 


- iK p(0) P^(r,0) 


= 0 


(r = 1) 


( 10 ) 


where is the two-dimensional Laplace operator in r and 0 and the eigenvalue 

\ is related to the axial wave number by the equation 



( 11 ) 


Equations (9) to (11) constitute an eigenvalue problem which must be solved to obtain 
the eigenfunction P«(r, 0) and the axial wave niimber An exact analytical solu- 

tion to this eigenvalue problem is possible and is developed in appendix A. The 
analytical solution requires finding the roots of a transcendental equation which 
include Bessel functions of many orders. To avoid the complexity of iteratively 
solving this transcendental equation, a numerical solution to this eigenvalue problem 
is obtained. The methodology for obtaining this numerical solution is described in 
the following paragraphs. 

The numerical technique for obtaining solutions to equations (9) and (10) is 
exactly as presented in reference 9 and only enough will be repeated to provide the 
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necessary nomenclature and ground work for the optimization method. Thus, the 
solution for the eigenfunction P^(r,0) is expressed in the form 


M 


N 


Pj^(r,e) = 


-EEO 

m=0 n=0 \ 


m,n 


cos m0 + B 


m,n 


sin m0 


J (X r) 
I m m,n 

' N 

m,n 


( 12 ) 


where ^ represents the hard-wall eigenvalues which satisfy the transcendental 

equation ^ normalization constant 

= r f cos\e J^(\ r) r dr d0 
nirn m m,n 

Galerkin's method allows determination of the modal coefficients ^ and ^ 

m,n m,n 

as well as the eigenvalue . 

There remains the problem of determining the number of circumferential and 
radial modes (M and N) which are needed in equation (12) to obtain an accurate 
resolution for the eigenmode Pp(r,0), Strictly speaking, N should be deter- 
mined such that the exact radial dependence of the exact eigenmode is accurately 
resolved by the series of hard-wall modes, or 

N 

J d J r) 

m X n m m,n 

n=0 ' 

where d^^ is an appropriately chosen constant. Studies reported in reference 9, as 
well as additional ones performed by the author, show that N = 9 gives good agree- 
ment between the exact and the numerically computed mode shape for a range of fre- 
quencies and m numbers considered here. Thus, results in this paper will be 
computed with N = 9. 

A simple rule of thumb for determining the upper index of summation M has been 
developed. This index should be chosen so that the incident wave can be accurately 
resolved by substituting z = 0 into equation (5), In addition, M is intimately 
related to the periodicity of the liner. For purposes of this paper, consideration 
is restricted to six circumferential modes. This does not mean that M = 5, since 
it is shown in appendix A that only selected values of m should appear in equa- 
tion (12) for any spatially discrete given source. Thus, only the first six of these 
selected values of m are included. This particular choice of M was based on con- 
vergence studies in which it was observed that virtually no change in the transmis- 
sion loss occurs as M is increased beyond this value for a range of frequencies and 
duct lengths. 


Scattering From Two Axial Discontinuities 


A finite-length circumferentially segmented liner combined with two semi- 
infinite hard -wall sections as depicted in figures 4 and 5 is now used for the study 
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of the problem of scattering from two axial discontinuities. An acoustic wave 
from the left is incident at z = 0, This results in a reflected wave in the 

first hard-wall section as well as a transmitted wave p^ in the second hard-wall 
section. It is necessary to compute the resulting transmission loss over the circum- 
ferentially segmented liner. In general, interaction of the incident wave with the 
peripheral liner is expected to result in modal transfer of acoustic energy among 
radial and circumferential modes of different orders. 

The constant geometry of the present problem makes it amenable to the mode match- 
ing technique (ref. 3). Thus, the incident, reflected, and transmitted waves in the 
hard-wall sections are expressed in terms of hard-wall duct modes 


M N 


P^(r,0,z) - Tv y. (aT cos m0 + b sin m0)-^^5 — — — — 


J (X r) exp(iK z) 
~ m,n 


m=0 n=0 


m,n 


(13) 


M N 


pN 


r,0,z) — ( 


m=0 n=0 


^ J (\ r) exp(-iK z) 

a^ cos m0 + b sin m0)-^!5^ — — — - 

m,n m,n 


(14) 


M N 


P^(r,0,z) = 7^ cos m0 + b^ sin m0} 

^ ^ m,n m,n ^ 


j ( \ r) exp[iK (z - L) ] 
m m,n m,n 


m=0 n=0 


N 


m,n 


(15) 


and 


2 2 2 
K = K - X 
m,n m,n 


(16) 


whereas the solution for the circumferentially segmented liner is expressed in terms 
of the circumferential liner duct modes 


P(r,e,z) = ^ {a exp(iKjj^z) + B exp[-iK^(z - L) ] } Pj^(r,0) (17) 

A=0 


in which the eigenmode Pj^ and the axial propagation constant are solutions to 

equations (9) to (11). The coefficients in the expansion for the incident wave are 
assumed to be known whereas those in the expansion for the waves in equations (14) 
to (17) have to be determined. It is instructive to note that would also be 

zero if the circumferentially lined section was infinitely long. (See eq. (5).) 
However, the finite termination of the lined section leads to reflections of the 
acoustic waves at z = L as indicated by the additional term exp[-iKj^(z - L) ] 
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in equation (17). Further, the length L has been included in equations (15) 
and (17) to assure that the waves running to the left and to the right attenuate 
properly in their direction of travel. 

Physical considerations require that the acoustic pressure and the corresponding 
axial derivative remain continuous at each junction between the circumferential liner 
and the hard-wall duct. Such requirements ensure continuity to the first order of 
mass and momentum between adjacent duct segments. Expressed here in vector form. 


{E} = {0} 


where {0} is a 4 x 1 null vector and {e} is a 4 x 1 residual vector, as shown in 
the following: 


^ 0 ^ 


{ 0 } = ( 




and 



where 

= P^(r,e,0) + p^(r,e,0) - p(r,0,o) 

I R 

Eg = P(r,e,L) - P^(r,e,L) 
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and 


=4 -■S""'®''’’ 


The mode matching equations are obtained by requiring that the residual vector be 
orthogonal to each of the functions Pj^(r,0), or 



{E(r,0)} P^(r,0) r dr d0 = {0} 


(Jl — Of If 2f •••f N) 


(18) 


where 


N=(2M+1)(N+1)-1 


Equation (18) constitutes a system of linear equations which can be expressed in 
matrix form as 


[G] {(J)} = {F} (19) 

where {(j)} contains the unknown coefficients of equations (14) f (15) f and (17) and 
{F} contains the known coefficients in the expansion for the incident wave, (See 
eq, (13),) Furthermoref the order of this system of equations is 4( 2M + 1 ) ( n + 1 ) , 
Equation (19) has been solved using Gaussian elimination to obtain {(()}, 

The transmission-loss function TL is still given by equations (7) and (8) for 
the finite circumferentially segmented liner as 


TL = 10 log4^Q[W(0)/W(L)] 


Howeverf the power of the incident wave is 


W(0) 



( 20 ) 


Equation (20) is a general expression which assumes that the incident acoustic wave 
is composed of multiple circumferential modes. Transmission-loss values presented in 
the "Results and Discussion” section of this paper are for incident waves consisting 
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of only a single circumferential mode, m = so that the power of the incident 

wave reduces to 


W(0) - (0) 

mg 


and 



(21 ) 


( 22 ) 


In addition, the output acoustic power is calculated from the transmitted acoustic 
wave for z = L as 


M 

W(L) = /^ W (L) 

m=0 


(23) 


and 


W (L) 
m 




(24) 


Now cutoff modes have the property Re(Kj^ ^) = 0 (i.e., K < ^m,n) whereas cut-on 

modes have the property >0* Tl^us, no acoustic energy is transported by 

cutoff modes. Hence, the summation in equations (20) to (24) actually reduces to 
summation for a finite number of cut-on modes. The transmission-loss function for a 
single circumferential -mode input becomes 


TL = 10 log 


10 


'M 

ni 

0 / m=0 


(25) 


whereas the transmission-loss function for circumferential mode m is 


TL 

m 


10 



(26) 


9 



Thus, TLj^ represents the amount of incident acoustic energy that is scattered and 
then absorbed from circumferential mode mg to circumferential mode m. It should 
be noted that no scattering between modes of different circumferential -mode orders 
occurs for a uniform or axially segmented liner. 


Optimization Method 

An optimization method was used to design and study circumferentially segmented 
liners for specified modal inputs and source frequencies. Although a variety of 
optimization methods exists, good success has been obtained with the so-called Quasi- 
Newton optimization algorithms. These algorithms combine the best features of 
steepest ascent and Newton's method. Apparently, the Davidon-Fletcher-Powell (DFP) 
optimization algorithm has been the most commonly used in duct acoustics (refs. 2 
and 3)* Thus, the DFP optimization algorithm with the transmission-loss function 
TL( the optimization parameter is used here. The free variables in the 
algorithm were taken to be the resistance R and the reactance ^ of the two 
admi ttance segments 


+ iy- 

1 

and 

6 ^ - 

P2 + iX2 


However, contour maps were first obtained to provide good starting values for the 
algorithm. More detailed information on the optimization algorithm used herein is 
given in appendix B. 

RESULTS AND DISCUSSION 

Trend studies for both uniform and circumferentially segmented liners are pre- 
sented. First, optimized segmented liners are presented and compared to optimized 
uniform liners for the semi -inf ini te duct depicted in figure 2. In order to reduce 
computational cost, optimization studies were performed only with L = 2. Second, 
studies are presented in which neither the uniform nor the segmented liner is opti- 
mized (i.e., off-optimum studies) with L = 2. These off-optimum trend studies 
include the effects of the two hard-wall sections (see fig. 4), so equations (23) 
and (24) are used to determine transmission losses. Finally, a cursory investigation 
of the effects of the lining length L is presented. These effects are investigated 
for the treated section installed within a hard-wall duct only. (See eqs. (23) 
and (24).) 
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Optimization Studies 


Optimization results for the semi -inf inite duct are presented in the following 
subsections. First, results for a plane wave input into a uniform liner are com- 
pared to^ results in an existing work to obtain confidence in the optimization 
method. Second, optimized segmented liners are compared with optimized uniform 
liners for a range of frequencies and source inputs. Segmented-liner studies 
consist of both hard-wall/soft-wall and soft-wall/soft-wall lining configurations. 

Plane wave incident into uniform liner. - Optimum transmission-loss values for a 
plane wave source are compared to those computed from the exact analysis of refer- 
ence 2 in figure 6 . Excellent agreement between the two analyses is observed. 

Optimum liner properties computed from the two analyses are presented in figures 7 
and 8 . Although the optimum resistance values compare favorably, agreement for the 
optimum reactance deteriorates with increasing frequency. This analysis shows 
excellent agreement at low frequencies but yields higher values than reference 2 as 
the frequency increases. However, the differences in reactance fail to change the 
optimum transmission losses. Thus, there may be more than one value of admittance in 
the complex admittance plane which gives the same transmission-loss value. 

Hard-wall/soft-wall optimization studies .- Figure 9 shows the optimum 
transmission-loss spectrum for a circumferentially segmented liner with different 
periodicities. The source is taken to be a plane wave with the optimum uniform liner 
used as a baseline for comparison. The width of the circumferential strips plays an 
important role in the liner performance. The performance of the segmented liner is 
an increasing function of the periodicity at each value of frequency. In addition, 
the uniform liner gives more suppression than the segmented liners. Note that the 
segmented-liner results approach the uniform-liner results as the periodicity 
increases. Optimum liner properties for the segmented liner are not presented. 
However, it was observed that for the segmented liners with high periodicity 

simply approaches the optimum admittance for the uniform liner. 

Optimxam transmission losses are presented in figure 10 for a segmented and a 
uniform liner. The incident wave for both cases was chosen to have a circumferential 
mode of m = 5 (Q(r, 0 ) = cos 50 ^ 5(^5 segmented liner has a periodicity 

of 5 so that the m numbers included in equation ( 12 ) are m=5, 10, 15, ... from 
the discussion in appendix A. As shown in figure 10, the uniform liner gives better 
performance than the segmented liner at the optimum point for all frequencies con- 
sidered. Optimum liner properties for the liners are shown in figures 11 and 12. 
Figure 11 shows an interesting trend for the optimum conductance. (The conductance 
is the real part of the complex admittance.) Note that the optimum conductance of 
the uniform and segmented liners are approximately equal. On the other hand, fig- 
ure 12 shows that the optimum susceptance of the two liners does not show such trends 
(i.e., the susceptance is the negative value of the imaginary part of the complex 
admittance) . 

We now turn to optimization studies at the higher frequencies for which high 
m-number sources are possible. A cursory investigation of the matrices of equa- 
tion (A9) reveals that the off-diagonal terms of these matrices increase in magnitude 
with increasing frequency, so that stronger coupling between the m numbers might be 
expected at the higher frequencies than the lower ones. Figure 13 compares the 
transmission losses for a uniform and a segmented liner. The source for both liners 
was a circumferential mode of m = 10 (Q(r,0) = cos 100 J^q^^IO 0^^^ with the 

periodicity of the segmented liner equal to 10 (T = 10). The segmented liner gives 
greater broadband performance than the uniform liner. Further, the performance 
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of the segmented liner approaches that of the uniform liner as the frequency gets 
closer to the cut-on value of the m number of the source. The optimum conductance 
and the optimum susceptance for the uniform and segmented liners are depicted in 
figures 14 and 15. Optimum conductance values for both the uniform and the segmented 
liner decrease in magnitude with increasing frequency; however, the segmented-liner 
values are higher than the uniform-liner values. Optimum susceptance values for the 
uniform liner approach zero at most frequencies. In contrast, the optimum suscept- 
ance values for the segmented liner vary with frequency from 0.29 to 0.62. 

Figure 16 shows a comparison of the optimum transmission-loss spectrums for a 
uniform and a segmented liner in which the circumferential mode order of the source 
is 15 (Q(r,0) = cos 150 Note that the segmented liner gives higher 

transmission losses than the uniform liner. Optimum conductance values for the 
uniform and the segmented liner are shown in figure 17. Results in this figure show 
that the conductance of the segmented liner is roughly twice that of the uniform 
liner. Optimum susceptance values for the two liners are given in figure 18 and show 
trends consistent with figure 13* Note that the optimum susceptance of the uniform 
liner is approximately zero, whereas that of the segmented liner is not. 

Soft-wall/soft-wall optimizati on studies.- The results presented so far 
restricted one section of the segmented liner to be a hard wall so that 
zero. Studies were also performed in which and were varied in the opti- 

mization algorithm. This lining configuration for which ^2 is not restricted to 
zero is referred to as the soft-wall/soft-wall circumferentially segmented liner. 
Further, results for the soft-wall/soft-wall segmented liner were computed for the 
same sources as for the hard-wall/soft-wall optimization studies. Optimum properties 
for the plane-wave source and the source with circumferential dependence of m = 5 
converge with the uniform-liner optimum properties where ” ^2* Although at some 

frequencies the admittances were different, the optimum transmission loss for the 
segmented liner was the same as the uniform liner to within 1 dB. A significant 
result was obtained for the sources with m = 10 and m = 15. For these sources, 
the optimization algorithm shows that p 2 = 0 + Oi was the optimum P 2 with p^ 
being identical to the value obtained from the hard-wall/soft-wall optimization 
studies. Further, for the frequencies in figure 13 for which the uniform liner gave 
slightly higher transmission losses than the segmented liner, the fully treated 
segmented-liner results converged with the uniform-liner values where p^ ” p2* 

Thus, a uniform liner stripped with aluminum tape may provide the best method of 
increasing the performance of a uniform liner for high m-n\imber sources. Further, 
removing the strips of aluminum tape and installing acoustic treatment may actually 
decrease the performance of the liner. 

Off -Optimum Studies 

In order to verify the accuracy of the mode matching equations used herein, 
sample calculations for a uniform liner with L = 2 were computed and compared to 
results obtained from the analysis of Zorumski (ref. 1). Computations obtained from 
Zorumski's analysis were based on the assumption that the superposition of 10 soft- 
wall modes could adequately represent the incident wave. In addition, the version of 
Zorumski *s analysis used included only single interface effects. Therefore, for 
comparison purposes, mode matching equations were developed for a single interface 
so that the second hard-wall section (see fig. 5) was removed and replaced by a 
nonreflecting termination at z = L. 
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Values from sample calculations are presented in figure 19 for an incident plane 
wave (P^(r,9, 0) =1). The admittance values used at each frequency to compute the 
transmission losses are optimum for no reflections at z = L. As depicted in fig- 
ure 19, excellent agreement between this analysis and that of reference 1 is obtained, 
A plot of values from a second sample calculation is presented in figure 20, Although 
the admittance values used to obtain the transmission losses in figure 20 are the same 
as those used in figure 19, the incident wave has been changed to 

P^(r,e,0) = cos 0 J (\ r) 

1 1 f u 


Discrepancies of about 11 percent at K = 4 and 26 percent at K = 5 are observed 
in the transmission losses computed from this analysis and that of Zorumski, Excel- 
lent agreement is obtained for K > 5. Further, it should be noted that by increas- 
ing N from 9 to 20, the discrepancy at K = 4 was reduced to less than 2 percent, 
though no significant improvement in the discrepancy at K = 5 was observed. It is 
believed this discrepancy can be removed by a superposition of more than 10 soft-wall 
modes to represent the incident wave in Zorumski’s analysis. Overall, good agreement 
has been obtained between this analysis and the analysis used in reference 1 . This 
leads to a degree of confidence in the mode matching equations developed in this 
paper. 


Having developed confidence in the analysis, we can now proceed to show that 
the analysis gives results consistent with the experimental observations of Mani 
(ref, 5). Figure 21 compares transmission-loss values for a uniform and a segmented 
liner but for an incident wave with m = 1 0 (P^(r,0,O) = cos 100 0^^^* 

values were computed for L = 2 with equation (25) being used to compute^ the 
transmission-loss spectrum. The uniform-liner admittance was set equal to , where 
the conductance and susceptance vary with frequency as shown in figures 13 and 14, 

The segmented-liner admittance function was obtained by setting T = 10, ^2 = ^ 

and the same as the uniform liner (i.e., the segmented liner is simply the 

uniform liner stripped with 10 evenly spaced strips of aluminum tape). As depicted 
in figure 21 , significant increases in performance can be realized by stripping the 
uniform liner, particularly at the frequencies closest to cut on of the m = 10 
circumferential mode. It must be emphasized that the increased suppression of the 
segmented liner over the uniform liner as depicted in figure 21 may not hold 
generally since the admittance p^ was deliberately chosen so that the segmented 
liner was optimized. Thus, different values of the admittance p^ may show the 
uniform liner to be superior to the segmented one. This may, however, be repre- 
sentative of the effect observed by Mani in reference 5, 

A more rigorous comparison of the performance of a uniform and a segmented liner 
is obtained when both liners are optimized, or tuned, at a fixed frequency. This 
allows a comparison of the off-optimum performance of the two liners. Figure 22 
compares the transmission-loss spectrum for a uniform and a segmented liner with a 
hard-wall/soft-wall admittance variation. Both the uniform and the segmented liner 
are tuned at K = 20 so that they are optimum only at this frequency. The segmented 
liner gives better performance not only at the tuning frequency, K = 20, but at 
other frequencies as well. Similar trends were observed when the tuning frequency of 
the two liners was K = 18 and K = 16* 
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Effects of Duct Length 


Figures 23 , 24, and 25 compare the transmission-loss spectrums computed from 
equation (25) for a uniform and a segmented liner for L = 1/2, 1, and 2. The 
incident wave for the figures was chosen to be P^(r,9,0) = cos 50 J 5 (X 5 q ^) 
whereas the admittance function for the segmented liner was chosen so that 
T = 5, p2 ^ ^1 “ 0,62i (i.e,, is an arbitrarily chosen 

value) . The admittance of the uniform liner was the same as for the segmented 

liner. Thus, the segmented-liner configuration can be obtained by stripping the 
uniform liner with tape. (This is the configuration proposed by Mani in ref. 5.) 

One can see from the figures that the segmented-liner performance increases relative 
to the uniform liner as dimensionless duct length increases. For L = 2, the seg- 
mented liner performs better than the uniform liner at all frequencies considered. 
This is surprising since only 50 percent of the segmented liner is treated (i.e., 

H ‘ »'• 

The total and modal transmission losses (see eqs. (25) and (26)) are probably 
also strongly affected by the duct length L. Figure 26 illustrates total and modal 
transmission losses obtained for a segmented liner as a function of dimensionless 
duct length. The computer calculations were made with K = 4.4, P 2 = 0 + Oi, 

P^ = 1 + 2i, T = 1, and P^(r,0,O) = cos 0 ( \^ 0^^* Equations (25) and (26) were 

used to compute the transmission losses so that tl^e effects of both hard-wall inter- 
faces are included. Results in the figure show that the circumferentially segmented 
liner has scattered acoustic energy into higher order circumferential modes than 
those present at the source. Since TL^ is at least 10 dB above TL 2 and TL^ , 
only a small portion of the total power is carried by the m = 3 mode in the second 
hard-wall section. Note also that TL^ and TL 2 are approximately equal when 
1 < L < 2. Thus, the circumferentially segmented liner has distributed energy 
equally between the m = 1 and m = 2 modes in the second hard-wall section for 
this range of L values. 


CONCLUSIONS 

Optimum and off-optimum properties of circumferentially segmented duct liners 
have been compared with those of uniform liners to analyze possible benefits of 
circumferential segmentation. Based on the results of this work, the following 
conclusions have been made: 

1 . Circumferentially segmented liners scatter acoustic energy among various 

circumferential wave numbers m. Results indicate that for some 
frequencies and duct lengths, the segmented liner may scatter source energy 
equally between a lower and higher order circumferential wave number. 

2. For low m-number sources, there is no advantage to an optimized circum- 

ferentially segmented liner over an optimized uniform liner. 

3. For high m-number sources, an optimized circumferentially segmented liner 

gives better performance than an optimized uniform liner. 

4. Segmented-liner configurations consisting of a hard-wall/soft-wall 

combination represent an optimum configuration for high m-number sources. 
Thus, replacing the hard-wall strips with acoustic treatment may decrease 
the performance of the liner. 


14 


5 


• Overall, the greatest benefit of circumferentially segmented liners over 
uniform liners occurred away from the optimum point. Off the optimum 
point, segmented liners give more effective broadband performance with a 
50-percent reduction in the amount of acoustic treatment. 

6. Increases in performance of the segmented liner relative to the uniform 

liner are greatest at frequencies closest to cut on of the circumferential 
mode number of the source. 

It should be emphasized that optimization studies were for a duct length-to~ 
radius ratio of 2, so these conclusions may not hold for all duct lengths. Results 
presented here are sufficiently encouraging (particularly for the high m-number 
sources) to warrant further studies involving circumferentially segmented liners in 
series or in checkerboard liners, the effects of multiple m-number sources, and the 
convective and refractive effect of a mean flow. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 2, 1982 
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APPENDIX A 


EXACT SOLUTION FOR CIRCUMFERENTIALLY SEGMENTED LINER 

This appendix presents the derivation of the exact solution to the elliptic 
partial differential equation (eq. (9)) 

and boundary condition 


Rj^(0) = 0 


(A1 ) 


where 


5P.(1 ,9) 

Rj^( 0) = — iK p(9) Pjj^d /0) 


(A2) 


and 


p(0) = 


P, + Pj 2(P, -Pj), 




sin T0 + 


sin 3T0 sin 5T0 




(A3) 


The goal of the derivation of this equation is to show how the periodicity of the 
acoustic liner T selects preferred values of circumferential mode orders. This is 
of great utility in that it provides greater physical insight into how source energy 
is redistributed into higher circumferential mode orders and it may lead to more 
realistic design criteria. 

It is easily, verifiable that each of the functions cos m0 J (X^r) and 
■f . m Jc 

sin m0 ^ solution to equation (9)* Thus, the most general solution to 

this equation is obtained by using the superposition principle 


P^(r,0) 




cos m0 + B sin m0 
m 


J [^0 


(A4) 


Although arbitrary values of the parameters A^, and satisfy equation (9), 

only selected values of these parameters will satisfy the boundary conditions imposed 
by the circumferentially segmented liner (eq, (A1)), Thus, it is still a matter of 
considerable mathematical difficulty to find the discrete values of the coefficients 
A^ and B^ and of eigenvalue which cause the boundary residual Rj^(0) to 

vanish as indicated by equation (A1 ) • 
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To satisfy the boundary condition imposed by the circumferentially segmented 
liner, one can make use of the fact that the set of functions (cos m0, sin m0) 
represents a complete set on the interval 0 to 2 it; and that only the zero function can 
be orthogonal to each member of a complete set* It is necessary that 


and 


r2% 

Jo 

r2it; 

Jo 


cos m0 d0 = 0 


sin m0 d0 = 0 


(m = 0, 1 , • • • ) 


(m — 1, 2, •••) 


(A5) 


(A6) 


Equations (A5) and (A6) are an alternate way of satisfying equation (A1 ) and are 
based on Galerkin's method (ref. 10). Further, the collection of equations (A5) 
and (A6) for each circumferential mode order m forms a set of simultaneous 
equations that can be expressed as 


CD] {X} = {0} 


(A7) 


where [D] is the complex coefficient matrix, {X} is the unknown vector of Fourier 
coefficients, and {0} is the null vector. Ordinarily, equation (A7) is solved to 
obtain the unknown vector {x} and the eigenvalue is determined by requiring 

that the determinant of the coefficient matrix vanish. It will now be shown that 
only a subset of this more general equation need be solved for any given source. 

The detailed structure of the matrix equation given by equation (A7) for a liner 
with periodicity T is 


[D] 


[Dq] 


[D.] 


T 



where [D^ ] is a square matrix, {Xg } is an unknown vector of the same order as 

[Dg], subscript s is the index 0, 1, •••/ T, and the superscript T denotes the 
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transpose. The diagonal character of [D] implies the periodicity of the acoustic 
liner decouples equation (A7) into (T + 1) eigenvalue problems 

= {0} 

[di m } = {0} 

[d2]{x2}={0> 1 

/ (A8) 


[d^]{x^} = {0} J 


where the vectors 


{Xg} are 


I A f A f ^ Am f •••/ 

•- s 2T-S 4T-S T-s ' 


T+s ' 3T-S ' 


] (s = 1, 2, T-1 ) 




/ [Af, A^„, sf, ..,1 

( L 0 2T 4T T 3T 5T -• 

b;5„, b^'„, b5„, ...] 

L 3rp' 5ff' ' 2T 4T 6T -* 


(s = 0) 


(s = T) 


where superscript Ji indicates the exact eigenfunction index. Thus, each sub- 
eigenvalue problem solved independently of each other and 

the eigenfunctions PpTr,0) for each subproblem are mutually orthogonal. This means 
that the only subproblems which need to be solved for any given source are those 
which contain the circumferential mode orders of the source. For example, only the 
subeigenvalue problem needs to be solved for an axisymmetric (cir- 

cumferentially uniform) source. Further, the matrices [^g] should be truncated at 
some finite order, (For example, results in the text use only the first six 
circumferential mode orders so that tD^] will be of order six,) Eigenvalues for 
each subproblem are obtained by setting each determinate in equation (A8) to zero 


det[Dg] =0 (s = 0, 1, T) 

so that these determinants determined the exact transcendental equations for a cir- 
cumferentially segmented liner. Note also that the exact transcendental equation for 
a uniform liner is obtained by setting allowing the periodicity to 

approach infinity. Under these conditions 

det[D ] = j’(Xo) - iKp. J (XJ (s = 0, 1, ,,,) 

S S X I S 

so that results here do degenerate to the exact solution for uniform liners. 
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DAVIDON-FLETCHER-POWELL (DFP) OPTIMIZATION ALGORITHM 


In this appendix, more detailed information concerning the optimization of the 
transmission-loss function is presented. This transmission-loss function is 
expressed in terms of an admittance vector TLCp^fp^) = TL{d}, where {d} is the 
4x1 admittance vector 



Starting with an initial admittance vector in which the superscript k 

is an iteration counter, it is the purpose of the DFP algorithm to move in a 
direction specified by the ascent vector so that a greater value of the 

transmission loss is obtained at the next value of the admittance vector } . 

The steps in the DFP algorithm (ref, 11) are as follows: 


1, Start with an initial admittance vector and any 4x4 symmetric 

positive-definite matrix (Here, was chosen as the identity 

matrix, ) 

2, Beginning with k = 0, compute where is the 

gradient vector of TL at 

3, Obtain the parameter a which maximizes 

TL({d^} + a*' {6^}) 


4, Set 

= {d’"} + a’" {6’"} 

{y’'} = 


and 


] . [H^] . Hf r - 




(y IK](y I 


where the superscript T denotes transpose, 
5, Set k = k + 1 and return to step. 1, 
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For this work, the iteration process was continued until the difference between 
l^k+1 I was at least 10~^. Note that the square matrix [h^] is 

updated at each iteration in an attempt to have the algorithm converge more rapidly 
than the method of steepest ascent (in which case [h^] is the identity matrix) and 
in an attempt to be computationally more efficient than Newton's method (in which 
case [h^] is the inverse of the matrix of second partials, or the inverse Hessian 
matrix, for which evaluation is either too costly or impractical) . 

Although the DFP algorithm has worked well in the past on duct transmission 
problems, there are still defects in the method which should be considered when 
optimizing acoustic liners. For example, steps 2 and 4 require the computation of 
the gradient of the transmission-loss function, TL({d^}). This gradient must be 
computed numerically and the error in the nximerical computation affects the 
convergence of the algorithm. In this paper the gradient was computed by finite 
differencing. On the other hand, step 3 of the DFP algorithm requires that the 
function TL({d^} + a^{6^}) be maximized to determine the search parameter 
However, since this function cannot be maximized analytically, numerical techniques 
must be used, which in turn leads to inaccuracies in the search parameter a . Such 
inaccuracies in can cause the algorithm to diverge or wander about meaningless ly 

in the admittance plane. The method of repeated parabolic interpolation was used in 
this paper to determine in step 3 of the DFP algorithm. Generally, it has been 

the author's experience that no major problems with the algorithm occur with the 
method of computing the parameter and transmission-loss gradient used 

here. A major problem did occur with the initialization vector because the 

DFP algorithm is not globally convergent. Originally, the optimum uniform values 
were used to initialize the algorithm; however, these values proved to be a local 
optimum, and larger transmission-loss values could be obtained by initializing the 
algorithm at nonuniform values. 

To alleviate the problems associated with the initialization of the DFP 
algorithm, an alternate method (henceforth referred to as the DPPC method) has been 
used. in the DFPC method, a locus of two-dimensional contours is used to obtain 
initialization values which are close to global optimums. These initialization 
values are then used as inputs to the DFP optimization algorithm. The steps in the 
DFPC method are the following: 

1. Vary R 2 and 1 to 10 in increments of 0.1. 

2. Obtain the optimum and transmission-loss value at each increment on 

or X 2 contour mapping. 

3. Obtain the global optimum from the loci of optimums obtained at each value 

of 

4. Now use the DFP algorithm with the global optimums obtained from step 3 as 

starting values. 

Although the DFPC method gives global optimum values, it will be quite expensive 
because a contour map must be obtained at each value of ^ 2 * Generally, results 
presented herein use 10 points on both the real and the imaginary axis of to 

obtain the contour maps in which the increment spacing in both directions was 
unity. However, at some frequencies the increment spacing had to be reduced further 
to obtain smooth contours, particularly at frequencies near cut on. 
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SYMBOLS 




amplitude of right and left moving wave in circumferentially 
segmented liner 




coefficients in exact expansion of eigenfunction P^(r,0) 





T5 Z 

®m,n 


a 


coefficients in hard-wall duct expansion of the exact eigen- 
function Pjj^(r,0) 

duct radius 


m,n' m,n 


m / n m / n 


a^ 

m,n' m,n 


c 


m 

{dM 

[D],[Dg],[G] 

{E} 


'® 2 '^ 3'^4 


{F} 

( 9 ^) 

[H^] 


i 


K 

L 


coefficients of Bessel-Fourier expansion of incident wave in 
first hard-wall section 

coefficients of Bessel-Fourier expansion of reflected wave in 
first hard -wall section 

coefficients of Bessel-Fourier expansion of transmitted wave in 
second hard-wall section 

ambient speed of sound 

coefficient in hard-wall duct expansion of Bessel function 
complex admittance vector 

complex admittance vector with iteration k 
complex coefficient matrices 
residual vector 

components of residual vector {E} 
known-constant vector 

gradient of transmission-loss function with iteration k 
symmetric positive -definite matrix of order k 
unit imaginary number 

Bessel function of first kind of order m 

wave number 

axial wave numbers 

predetermined length over which transmission loss is computed 
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M,N 


total n\imer of m and n values 


m 


— Of •••/ M— 1 


N _ normalization constant for hard-wall basis functions 

n 


P(r, e,z) 

P^(r, e,z) ,P^(r, 0,2) 

Pj^(r, 0) 

P^(r, 0,z) 

Q(r, 0) 

R 

Rb(0) 

Re( ) 

r,z 

T 

t 

W(z) 

W^(z) 

{xjriXg} 

p( 0) , , p2 

Pav 

{ 6 ’^} 

V,n 

p 

{<!.} 


= 0, 1, N - 1 

acoustic pressure field in duct with circumferentially segmented 
liner 

incident and reflected pressure field in first hard-wall 
section 

acoustic pressure eigenfunction 

transmitted wave in second hard-wall section 

source pressure function in semi -infinite duct 

resistance of segmented liner 

boundary residual 

real part of complex expression 

radial and axial coordinates 

periodicity of segmented liner 

total transmission-loss function 

transmission loss of circumferential mode m 

time 

total sound power at axial location z 

sound power in circumferential mode m at axial location z 
unknown-constant vectors 
search parameter 
acous tic admi ttances 

average admittance for circumferentially segmented liner 
ascent vector with iteration k 

eigenvalues of circumferentially segmented and hard-wall duct 
ambient density 
unknown vector 
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0 ) 

X 


angular frequency 
reactance of segmented liner 
Laplace operator 


Subscripts: 

Ji exact eigenfunction index 


m 


circumferential wave number index 


n 


radial eigenfunction index for hard-wall duct 


s subeigenvalue problem index^ s=0^ 1, T 

Superscripts: 

k iteration index 


T transpose operator 

A bar over a symbol denotes a dimensional quantity. A prime denotes derivative 
with respect to the argument. 
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First hard-wall 
section 


r 


Circumferentially 
segmented liner- 


Second hard-wall 
section -j 



Figure 4.- Circumferentially segmented liner installed in hard-wall duct 



Figure 5.- Configuration used to develop mode matching equations. 





Figure Optimum transmission-loss spectrum for plane wave incident 

into an infinite uniform liner. 
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Result from this analysis 
Result from ref. 2 
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Dimensionless frequency, K 


Figure 7.- Optimum resistance for transmission-loss spectrum of figure 6 







Figure 9.- Comparison of optimum transmission-loss spectrums for uniform 
and segmented liners for plane-wave source. 
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Figure 10«- Comparison of optimum transmission loss for standing-wave source 








Uniform liner 



I I I I I I I 

7 18 19 20 21 22 23 24 

Dimensionless frequency, K 

Figure 14.- Optimum conductance for standing-wave source, m = 10. 








Optimum susceptance 



Figure 17.- Optimum conductance for standing-wave source, m 
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Uniform liner 
Segmented liner, T = 1! 


- L _ r^ I - I I- : 1 

26 28 30 32 

Dimensionless frequency, K 


Figure 18.- Optimum susceptance for standing-wave source, m 






Dimensionless frequency, K 


Figure 20 •— Comparison of transmission losses predicted from two different 
analyses for standing— wave source, m = !• 






Dimensionless frequency, K 


Figure 22.- Off-optimum performance properties for uniform and segmented 

liners tuned at K = 20. 
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Dimensionless frequency, K 


Figure 23.- Transmission-loss spectrums for uniform and segmented (stripped) 
liners for standing-wave source. m = 1; L = 1/2. 




Dimensionless frequency, K 


Figure 24.- Transmission-loss spectrum for uniform and segmented (stripped) 
liners for standing-wave source. m = 1 ; L = 1 . 






Dimensionless duct length, L 


Figure 26 Total and modal transmission losses for segmented liner with 

standing- wave source* m = 1 . 
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